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Abstract 

The transport of high energy heavy (HZE) ions through bulk ma- 
terials is studied with energy dependence of the nuclear cross sections 
being neglected . A three-term perturbation expansion appears to be ad- 
equate for most practical applications for which penetration depths are 
less than 30 g/cm 2 of material. The differential energy flux is found for 
monoenergetic beams and for realistic ion beam spectral distributions . 
An approximate formalism is given to estimate higher order terms. 
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1. Introduction 

Although heavy ion transport codes for use in space applications are in a relatively advanced 
stage at Langley Research Center (refs. 1 through 3), it seems prudent to further develop the 
theory for comparison with laboratory experiments, which has been recently neglected since 
initial efforts began several years ago (refs. 4 through 7). In the present report, we begin with 
the most simplified assumptions for which the problem may be solved completely. Solutions to 
a more complete theory may then be compared with prior results as limiting cases. In this way, 
the more complete but approximate analysis will have some basis for evaluating the accuracy 
of the solution method. The lowest order approximation will be totally energy-independent. 
The most complicated solution to be considered herein will have energy-independent nuclear 
cross sections but will treat the energy-dependent atomic/molecular processes and the energy 
spread of the primary beam. A fully energy-dependent theory must await further development, 
although some terms have been previously evaluated (ref. 6). 

2. Energy-Independent Flux 

If the ion beam is of sufficiently high energy that the energy shift due to atomic/molecular 
collisions brings none of the particles to rest in the region of interest, then 


r d 1 


4>j(x) — ^ ^ m jk&k fiki 3 ') 

k 


( 2 . 1 ) 


where <f>j(x) is the flux of type j ions, a j is the nuclear absorption cross section, and is 
the fragmentation parameter for producing type j ions from type h. The solution for a given 
incident ion type J is given in terms of a set of ^-functions as follows: 


g(j l) = exp (-<T jx x) 


9 {j h 32) * * • •> Jn-i jn+l) 


9 {jl 2 32 1 • • • t Jn—l i 3 ti) 9 (jl •> • • * ? 3n— 1 1 Jn+1 ) 

&jn+l — &jn 


for which the solution for the type j ion flux is written as 


4>f\x) = 6j. g(j) 

,(1 ), , n exp(-Sjx) - exp(-Sjx) 

<t>) (x) = rrijjcrj g{j, J) = rrijjaj J - 

J vj ~ Cj 

4>f\x) = 53 m jk a k m k jcrj g(j, k, J) 

k 

<pf\x) = 53 m jk°k Wkivimijvj g(j,k,l,J) 

k,l 

with 

M x ) = ^( x ) 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


This solution is equivalent to that derived by Ganapol et al. in reference 8. We now consider 
some applications of the above formalism. The cross-section data base is discussed elsewhere 
(ref. 9). 

2.1. Neon Beam Transport 

We first note in the case of 20 Ne incident on water that 19 Ne and 19 F have only one 
contributing term in equation (2.8). These are shown in figure 2.1: Also shown in figure 2.1 



are the fluxes of various isotopes of secondary ion fragments. The effect of successive terms of 
equation (2.8) is shown in table 2.1 for the lS 0 flux. It is clear from the table that the fourth 
and higher order collision terms are completely negligible and that third collision terms are a 
rather minor contribution. Hence, a three-term expansion as we have used in the past (refs. 5, 
6, and 7) appears justified. The relative magnitude of the terms contributing to the 'Li flux 
generated by the 2 ^Ne beam is presented in table 2.2. The fourth collision term is negligible at 
small penetration distances and small, but not negligible, at distances greater than 30 cm. The 
greater penetrating power of the lighter mass fragments is demonstrated in figure 2.2. Also 
note the difference in solution character due to the importance of the higher order term. 

2.2. Iron Beam Transport 

We first note in the case of 56 Fe incident on water that 55 Fe and 55 Mn have only one 
contributing term in equation (2.8). The 54 Mn has but two terms, and the slight difference in 
solution character can be seen in figure 2.3. Results for 52 V are also shown. The convergence 
rate of equation (2.8) is demonstrated in table 2.3. Again we see the fourth collision term to be 
negligible, whereas the three-term expansion we have used before seems quite accurate at these 
depths for these ions. In distinction to prior results, the 16 0 flux has significant contributions 
from higher order terms for depths beyond 20 cm as seen in table 2.4. Clearly, a more complete 
theory using higher order terms is required than previously used for ion beams of particles 
heavier than 20 Ne. The different solution character of the lighter mass fragments is clearly 
demonstrated in figure 2.4. 


Table 2.1. Normalized Contributions to 15 0 Flux From Successive 
Collision Terms for 2 ^Ne Transport in Water 


Fragment 

term 

15 0 flux at x of — 

10 cm 

20 cm 

30 cm 

40 cm 

50 cm 


1.00E0 

1.00E0 

l.OOEO 

1.00E0 

l.OOEO 


1.01E-1 

2.01E-1 

3.02E-1 

4.03E-1 

5.04E-1 

0<3) 

2.63E-3 

1.05E-2 

2.36E-2 

4.18E-2 

6.52E-2 

4.W 

3.31E-5 

2.52E-4 

8.58E-4 

2.03E-3 

3.95E-3 


Table 2.2. Normalized Contributions to 7 Li Flux From Successive 
Collision Terms for 2 ®Ne Transport in Water 



7 Li flux at x of— 

Fragment 






term 

10 cm 

20 cm 

30 cm 

40 cm 

50 cm 

0<n 

l.OOEO 

l.OOEO 

l.OOEO 

l.OOEO 

l.OOEO 

0(2) 

1.62E-1 

3.20E-1 

4.72E-1 

6.18E-1 

7.58E-1 

0(3) 

1.15E-2 

4.53E-2 

9.98E-2 

1.73E-1 

2.63E-1 

0(0 

4.02E-4 

3.16E-3 

1.04E-2 

2.39E-2 

4.53E-2 
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Table 2.3. Normalized Contributions to 52 V Flux From Successive 
Collision Terms for 56 Fe Transport in Water 


Fragment 

term 

52 V flux at x of— 

10 cm 


30 cm 

40 cm 

50 cm 

— 0d) 

1.00E0 

1.00E0 

1.00E0 

1.00E0 

1.00E0 

0(2) 

7.91E-2 

1.52E-1 

2.37E-1 

3.15E-1 

3.94E-1 

0(3) 

2.37E-3 

9.48E-3 

2.13E-2 

3.79E-2 

5.91E-2 

0(4) 

2.24E-5 

1.73E-4 

5.93E-4 

1.41E-3 

2.75E-3 


Table 2.4. Normalized Contributions to 16 0 Flux From Successive 
Collision Terms for 56 Fe Transport in Water 


Fragment 

term 

16 0 flux at x of * 

10 cm 

20 cm 

30 cm 

40 cm 

50 cm 

0OT 

1.00E0 

1.00E0 

1.00E0 

1.00E0 

1.00E0 

0(2) 

5.87E-1 

1.12E0 

1.59E0 

2.00EO 

2.36E0 

0(3) 

1.86E-1 

7.08E-1 

1.49E0 

2.46E0 

3.56E0 

0(4) 

3.06E-2 

2.63E-1 

9.44E-1 

2.33E0 

4.72E0 



Figure 2.1. Ion fragment flux of various isotopes as a function of depth in water for a 20 Ne 
incident beam. 


3 














Figure 2.4. Flux of light ion fragments as a function of depth in water for a 56 Fe incident beam. 

3, Monoenergetic Ion Beams 

In moving through extended matter, heavy ions lose energy through interaction with atomic 
electrons along their trajectories. On occasion, they interact violently with nuclei of the matter 
and produce ion fragments moving in the forward direction and low energy fragments of the 
struck target nucleus. The transport equations for the short range target fragments can be 
solved in closed form in terms of collision density (refs. 5 and 6). Hence, the projectile 
fragment transport is the interesting unsolved problem. In previous work, the projectile ion 
fragments were treated as if all went straightforward (ref. 4). We continue with this assumption 
herein, noting that an extension of the beam fragmentation model to three dimensions is being 
developed (ref. 9). 

With the straightahead approximation and the target secondary fragments neglected (refs. 4, 
5, and 6), the transport equation may be written as 


'd_ 

dx 



(E) + CTj 


<f>j (x, E ) 


mjkVk <t>k( x >E) 

k>j 


( 3 - 1 ) 


where <f)j(x , E ) is the flux of ions of type j with atomic mass Aj at x moving along the x-axis at 
energy E in units of MeV /amu, aj is the corresponding macroscopic nuclear absorption cross 

section, Sj(E) is the change in E per unit distance, and rrijk is the fragmentation parameter 
for ion j produced in collision by ion k. The range of the ion is given as 



dE f 

Sj(E') 


(3.2) 


The stopping powers used herein are based on Ziegler’s fits to a large data base (refs. 11 
through 16). There is some controversy as to the stopping powers to be used (ref. 10). The 
analysis in reference 10 was biased to the stopping power used in the data analysis and is the 
same as that used in the PROPAGATE code. The values in the HZESEC computer code are 
similar to those in PROPAGATE. Neither the PROPAGATE nor HZESEC stopping powers 
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have been compared with the data base collected by Ziegler as far as is known to us. We 
continue to use Ziegler’s work until more definitive comparisons compel us to do otherwise. 

The solution to equation (3.1) is to be found subject to boundary specification at x = 0 and 
arbitrary E as 


<i> j (0,E) = F j (E) 

Usually Fj(E) is called the incident beam spectrum. 

It follows from Bethe’s theory that 


Sj(E) 




S P (E) 


for which 


Z i Z v 

-t M E ) = / ME) 


(3.3) 


(3.4) 


(3.5) 


The subscript p refers to proton. Equation (3.5) is quite accurate at high energy and only 
approximately true at low energy because of electron capture by the ion which effectively 
reduces its charge, higher order Born corrections to Bethe’s theory, and nuclear stopping at the 
lowest energies. Herein, the parameter vj is defined as 


Z? 


A 3 


so that 


Vj Rj{E) = v k R k (E) 


(3.6) 


(3.7) 


Equations (3.6) and (3.7) are used in the subsequent development, and the energy variation in 
vj is neglected. The inverse function of Rj(E) is defined as 


E = RJ 1 [RjiE)] 


(3.8) 


and plays a fundamental role subsequently. For the purpose of solving equation (3.1), define 
the coordinate transformation (refs. 5 and 6), 


7]j = x~ Rj{E) | 
= x Rj(E) j 


and new functions 


where 


= Sj{E) (t>j{x 
Xk(vj,Zj) = Xk(vk^k) 

Zj + Vj = Zk + *tk 

V] - Zj = ^(.Vk - Zk) 


,£) j 


(3.9) 


(3.10) 


(3.11) 


for which equation (3.1) becomes 

d 


+ ) Xjbi'Zj) = E m j^k-^ Xk(Vj,Zj) 


(3.12) 
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where the <jj is assumed to be energy independent. Solving equation (3.12) by using line 
integration with the integrating factor, 


VjiVjitj) = exp 



(3.13) 


results in 


Xj(rtj,Zj) = exp + Tjj ) 


Xj ( i €j ) 


l Hi 1 1 ] v 

+ o / , exp ~ V) E Xk^Kk) dr l (3-14) 


where 


,/ i = !2L+S,' + !2ULf2 { . 

‘ 2>'* 2i/ t '' 


d - V k~ v j _r , v k + Vj 

(i - — " + -sr ^ 


and the boundary condition (eq. (3.3)) is written as 

Consider a Neumann series for equation (3.14) for which the first term is 

X?\vj,Zj) = 


and the second term is 


(3.15) 


! v <j) = ex P + (j) Sj K _1 (Ci)] Fj [fi- 1 ^)] (3.16) 


xf ] (Vj,Sj) = \ e*p[|o>-(i/ - ty)] E m >fc ff fc^exp[-^fc(»?J: ~ £*) 

n[V(«*)] dT >' (3-17) 

An expression for Xj {Vj^j) is derived once equation (3.17) is reduced and higher order terms 
can be found by continued iteration of equation (3.14). These expressions (eqs. (3.16) and 
(3.17)) are now simplified for a monoenergetic beam of type M ions. 

The boundary condition is now taken as 


F j (E) = 6 jM 6(E-E 0 ) (3.18) 

where 6jj \f is the Kronecker delta, 6{ ) is the Dirac delta, and E 0 is the incident beam energy. 
Thus, 

Xj(-Zj,Sj) = Sj [VM S JM 6 [^7 1 (e j )--E 0 ] = 6ju S^j-RjiEo)} (3.19) 
for which becomes 


xf\vj,Zj) = 6j M ex p -^<7j(T)j+Sj) I Sfc-RjiEo)} 


(3.20) 
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and Xj becomes 


xf\vj,Zj) = \ m jM a M ^-ex p ~ + 


x 6 [^'m~ R m( e o)] dr]' 


(3.21) 


where is given by equation (3.15) for k = M. The contribution to the integral (eq. (3.21)) 
occurs at _ , „ 


VVV/UIO 0 . 

I'M - Vi I'M - ", 

provided that rf lies on the interval — < rf < r}j so that 




(3.22) 


(3.23) 


The simplified form in equation (3.23) may now be used to calculate the next iteration of 
equation (3.14): 

xfHvj^j) = \ £ rn jk a k m k M^M Vk j ex P + ft 


where 


- 2 a k (4 -*?) - toi - 4) *n" 


„ _t / k + Vj n Vk v j , 
Vk - 2u k V + 2u k 


/, _vk~ n . n+^j , 
** " 2u k V + 2v k 


(3.24) 


(3.25) 


UiiU . 

_ _2*m_ fijf(Eo) _ SIL + jStf (3.26) 

I'M - V k V M ~ u k 

with the requirement that -£ k <fj< v k - The inverse of the transformation is now applied to 
obtain from equation (3.20) 


V M ~ Vk 


(3.26) 


<j>f\x,E) = ~ j -- exp(-gjx) 6j M 6[x+Rj(E)-R M (E 0 )] 


Sj(E) 


(3.27) 


and from equation (3.23) 


so long as 


^ 1)(X ’ E) = I^E) miM<TM WM~- vj\ eXp {~h [X ~ Rj{E) ~ V ' ] 

-\ a M [ x + R j( E ) + »/]} 

— [Rm(Eo) - *] < Rj(E) < ^ R M (Eo ) - x 


(3.28) 


(3.29) 


Otherwise <j>^\x,E) is zero. After a complicated but straightforward manipulation, a similar 

^ ( 2 ) 
result may be obtained from equation (3.24) for <jr- (x,E). 
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In reducing equation (3.24), it is useful to define 


X M — ^ {ik + v) 

(3.30) 

x k=\ (v'k - v) 

(3.31) 

x i = \ (Vj - v") 

(3.32) 


and make a change in variables as 

,(2)/ 77. \ _ 1 Vj [ X 3 U 

<pj (x,&) - 2s m jk <r k m kM < ^ M j x exp(-cr M i M - a k x k - crjXj) dxj 

(3.33) 

where the integral is understood to be nonzero only in the physically allowed regions as 
explained presently. One may easily demonstrate 


X M + x k -|- Xj — X 

VM-^M d" V k X k d” VjXj = V M R-m{Eo) ~ Vk R-k{E) 

for which the parametric solution is given as 

_ _ V M Rm(Eo) ~ v k [R k (E) + x] + ( v k - Vj)xj 

X M — 

v M~Vk 


x _ ^m\Em{E) + x\--V}4 Rm{Eq) - (i>M ~ Vj)xj 

k v M -v k 

The requirement that x and x k be bounded by the interval 0 to x — xj yields 


(3.34) 

(3.35) 


(3.36) 

(3.37) 


° 1 


f X ' 

j v k [R-k(E) d- x] - V M R M (Eo) | 

> < Xj < i 

1 v M [R M {E) +x]~ v m R m (E 0 ) > 

l Vk - Vj J 


l v M - V j 


as the appropriate limits for the integral in equation (3.33) when > v k > Vj. In the braces 
in equation (3.38), we always choose the most restrictive value for the limit. The requirement 
of equation (3.38) also implies the result that 


RJf [R m (E 0 ) — x]<E<R jT 1 


V M Rm(Eq) - VjX 

v k 


(3.39) 


as the range over which the result of equation (3.33) is not zero. In the event that u k > v M > 1 
then 


° 1 f 

YM. tSEl t x ) ~ V M Rm(Eo) f - ^ I 'k [ Rk(E ’) + x]-v M Rm{Eq ) 

V M ~Vj ) l v k - Vj 

As a result of equation (3.40) 


(3.40) 




V M 

L ^ 


Em(E 0 ) ~ x 


<E<R~J 


Rm(E 0 ) — i- x 

V M . 


(3.41) 
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In the event that vm > v j > follows that 


0 < x j < < 


[Rm{ E ) + x ~ Rm( E o)] 

V M - Vj 

umRm[Eq) - Vk Rk{E ) - t/fcx 

^ - n 


(3.42) 


where the lesser of the three values in the braces is used as the upper limit of Xj for which the 
integral of equation (3.24) is not zero. As a result of equation (3.42) 


R a / [Rm(Eq) ®] ^ E — R-k 


- A -R M {Eo)-x 
v k 


(3.43) 


The integral in equation (3.33) may now be evaluated as 






Sj(E) \vm v k\^jkM 
exp( <tm x Mu &k x ku — &j x ju)\ 


[exp {-(Tm x MI - a k x kl - Vj x jl) 


(3.44) 


where x\f u , xj. u , xj and x^\ are the values of equations (3.36) and (3.37) evaluated at the 
corresponding upper and lower limits of Xj and 


AjfcM = a j + 


in ~ vj) 


(I'M ~ Vj) 


WvM-Vk)*™ (VM-Vk) k 


(3.45) 


Higher order terms are similarly derived. 

The total integral flux associated with each term may be evaluated as 


*«(*) = J™ 4>f(x,E) dE 


One may easily show that 

L 


00 AD, r*\ a jM [exp(-Cji) - exp(-<7 M x)] 

<P) ( x t E) aE = — — 

J (T\{ (Tj 


(3.46) 


(3.47) 


in agreement with equation (2.5). Furthermore 



°jk<nc M 

\ V M ~ v k\^jkM 



[exp(— <7jx) - exp(-<r A .x)] 


- [exp(-CTjx) - exp(-<r M x)]| (3.48) 

which agrees with equation (2.6) as — < > I'm- This relation of equation (3.48) and 

equation (2.6) has been used previously (ref. 7). 

3.1. Total Flux Comparisons 

The results of equations (3.28) and (3.44) are integrated numerically over their entire energy 
spectrum and given along with values from corresponding energy-independent solutions in 


table 3.1. The primary beam was taken as 20 Ne at 1380 MeV/amu. Clearly, the energy- 
dependent solutions appear quite accurate. 

3.2. Monoenergetic Beam Results 

The fluorine spectral flux seen at various depths in a water column is shown in figure 3.1. 
The primary beam was 20 Ne ions at 600 MeV/amu corresponding to a range of 30 cm. There 
is a clear structure due to the fluorine isotopes shown in the spectrum. The most energetic ions 
are 19 F. The 18 F and 17 F spectral components are clearly resolved. Only the 19 F is able to 
penetrate to the largest depth represented (35 cm). A similar, but more complicated, isotopic 
structure is seen in the oxygen spectra of figure 3.2. The greater number of oxygen isotopes 
contributing has a smoothing effect on the resultant spectrum. This effect is even more clearly 
seen in figure 3.3 for the nitrogen isotopes. Some of the smoothness results from the higher 
order term (p^ in the perturbation expansion. The boron flux of figure 3.4 shows very little 
isotopic structure. Qualitatively, similar results are obtained for an iron beam of the same 
range (30 cm) as shown in figures 3.5 to 3.9. 


Table 3.1. Total Flux From Energy-Independent Solution and Numerically Integrated Differential Spectrum 
[Values in parentheses are from energy- independent solution] 


Fragment 

Term 


Flux, cm 2 , at water depth x of- 

- 

5 cm 

20 cm 

18p 

00) 

0.00727 

(0.00717) 

0.01148 

(0.01140) 


4>( 2 ) 

0.00018 

(0.00018) 

0.00114 

(0.00114) 

17 0 

*<»> 

0.00729 

(0.00729) 

0.01173 

(0.01174) 



0.00017 

(0.00017) 

0.00112 

(0.00112) 

16q 

itn 

0.01350 

(0.01349) 

0.02193 

(0.02202) 


4>( 2 ) 

0.00029 

(0.00029) 

0.00191 

(0.00190) 

“N 

*(>) 

0.00470 

(0.00481) 

0.00796 

(0.00796) 


0< 2 ) 

0.00032 

(0.00033) 

0.00220 

(0.00220) 

13 C 

/ 1) 

0.00511 

(0.00521) 

0.00894 

(0.00887) 


<f>( 2 ) 

0.00032 

(0.00033) 

0.00224 

(0.00224) 

12 C 


0.00668 

(0.00682) 

0.01173 

(0.01178) 



0.00056 

(0.00056) 

0.00398 

(0.00398) 

1J B 

00 ) 

0.00417 

(0.00417) 

0.00735 

(0.00732) 


0(2) 

0.00036 

(0.00036) 

0.00259 

(0.00259) 
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Figure 3.1. Fluorine flux spectrum produced by a 20 Ne beam at 600 MeV/amu in a water 
column at various depths. 



Figure 3.2. Oxygen flux spectrum produced by a 20 Ne beam at 600 MeV/amu in a water 
column at various depths. 



Figure 3.3. Nitrogen flux spectrum produced by a 20 Ne beam at 600 MeV/amu in a water 
column at various depths. 



Figure 3.4. Boron flux spectrum produced by a 20 Ne beam at 600 MeV/amu in a water column 
at various depths. 
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E, MeV/amu 


Figure 3.5. Manganese flux spectrum produced by a 56 Fe beam at 1090 MeV/amu in a water 
column at various depths. 



Figure 3.6. Chromium flux spectrum produced by a 56 Fe beam at 1090 MeV/amu in a water 
column at various depths. 




Figure 3.7. Vanadium flux spectrum produced by a 56 Fe beam at 1090 MeV/amu in a water 
column at various depths. 



Figure 3.8. Fluorine flux spectrum produced by a 56 Fe beam at 1090 MeV/amu in a water 
column at various depths. 
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Figure 3.9. Carbon flux spectrum produced by a 56 Fe beam at 1090 MeV/amu in a water 
column at various depths. 


4. Realistic Ion Beams 

In the previous section, we assumed that a monoenergetic beam was present at the boundary. 
We now take the incident ion beam flux to be 





(E - £ 0 ) 2 

2A 2 


(4.1) 


where E 0 is the nominal beam energy and A is related to the half-width at half-maximum. 
The full solution is then found as a superposition of results from the previous section. The 
uncollided flux is found to be 


4> l l$(x,E) 


S m (E) 


exp {-cr M x) 



exp 


(Eg ~ Ernf 

2A 2 


(4.2) 


where R^(E m ) — R^{E) + x. One similarly arrives at 


<t>f\x,E) = ~ 


a jM v j 


Sj(E ) | vm 


— — exp |-^crj [x - Rj(E) - t]' 0 ] - [x + Rj(E) + j/ 0 ] 




where 




, _ _ VM+J^ J + T ] 

VM- Vj VM-Vj 


} 


(4.4) 

(4.5) 

(4.6) 
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The second collision contribution to the ion energy spectrum is similar: 


<t>f{x,E) = Y, 




where 


~k S j( E ) \ V M ~ v k\ &jkM 

M x Mu ^k x ku °j x ju)\ 

“(VHfii?)] 


[exp {<?M x Ml ~ a k x kl ~ cr j x jl ) 


1 

x - 
2 


Ei={ 


r m 


R 


M 


R M 


VkRkjE) + VjX 

»M 

\ Rm{E) + vjx ' 

V M 

v k (R k (E) + x) 


V M 


( U M > v k> Vj) 
(n > V M > Vj) 

(v M > t'j > v k ) 


(4.7) 


(4.8) 


R M [ r m(E) + x] 


E u =< 


R M 


\ l/ ki R k(E) + x) 


V M 

R \i [ r m(E) + x] 


{v M > Vk> v j) ) 
i v k > U M > Vj j) 

(V M > Vj > v k ) ) 


(4.9) 


and x M and x k evaluated at the upper and lower limit values of xj are obtained from 
equations (3.31) and (3.32). 

The elemental flux spectra were recalculated for 20 Ne ions at 600 MeV/amu with a 
0.2-percent energy spread assumed for the primary beam. The resulting fluorine flux is shown 
in figure 4.1. Although the spectral results are quite similar to the monoenergetic beam case, 
there is a considerable smoothing of the total spectrum. Similar results are obtained for the 
oxygen flux as well in figure 4.2. In distinction, the nitrogen and carbon spectra show only 
slight isotopic structure as seen in figures 4.3 and 4.4. Qualitatively similar results are obtained 
for the 56 Fe realistic beam as shown in figures 4.5 through 4.9. 
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E, MeV/amu 


Figure 4.1. Fluorine flux spectrum produced by a 20 Ne beam at 600 MeV/amu with a 
0. 2-percent energy spread in a water column at various depths. 



Figure 4.2. Oxygen flux spectrum produced by a 20 Ne beam at 600 MeV / amu with a 0.2-percent 
energy spread in a water column at various depths. 






E, MeV/amu 


Figure 4.4. Carbon flux spectrum produced by a 20 Ne beam at 600 MeV/amu with a 0.2-percent 
energy spread in a water column at various depths. 


E, MeV/amu 


Figure 4.5. Manganese flux spectrum produced by a 56 Fe beam at 1090 MeV/amu with a 
0. 2-percent energy spread in a water column at various depths. 



Figure 4.6. Chromium flux spectrum produced by a 56 Fe beam at 1090 MeV/amu with a 
0.2-percent energy spread in a water column at various depths. 




Figure 4.7. Vanadium flux spectrum produced by a 56 Fe beam at 1090 MeV/amu with a 
0.2-percent energy spread in a water column at various depths. 



Figure 4.8. Fluorine flux spectrum produced by a 56 Fe beam at 1090 MeV/amu with a 
0.2-percent energy spread in a water column at various depths. 
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Figure 4.9. Carbon flux spectrum produced by a 56 Fe beam at 1090 MeV/amu with a 
0. 2-percent energy spread in a water column at various depths. 

5. Approximate Spectral Solutions 

In the previous sections, the spectral solutions of the secondary ion flux were derived 
to second-order collision terms. Such a three-term expansion is not always an adequate 
representation of the transport solution. In this section, we derive approximate expressions 
for the perturbation series. Clearly, the more accurate results would be used to the order 
to which they are known, and the higher order terms would be taken to the approximate 
expressions of this section. 

5.1. Approximate Monoenergetic Beams 

The uncollided beam solution is taken as 

tf(x, E) = #(*) + flj(E) - *<(«.)] < 5 -D 

which is equal to the result in equation (3.27). The first-order collision term is approximated 
by noting that the energy dependence of the exponent of equation (3.28) is slowly varying in 
energy (refs. 4 through 7) resulting in 
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Similarly, 
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where Ej u and Eji are given by equations (5.6) and (5.7). In all the expressions for 0‘ n ' 1 given 
by equations (5.2), (5.5), and (5.8), the flux values are taken as zero unless 


Eji <E< Ej u (5.9) 

The approximate monoenergetic beam solutions are given in figures 5.1 through 5.4 and should 
be compared with the solutions found in section 3. The 17 0 flux at 20 cm of water is shown 
in figure 5.1 as contributed by the first collision term. The trapezoidal (solid) curve is the 
exact solution for the first collision term derived in section 3. The rectangular (dashed) curve 
is the approximate first collision term of equation (5.2). Terms for other fragment spectra are 
similar to those shown in figure 5.1. The solution for the second collision contribution to the 
O flux at 20 cm of water is shown in figure 5.2. The nearly rectangular solution (dashed 
curve) is the approximation given by equation (5.5). A triangular spectral function of the 
same energy interval could yield improved results. The spectra of fragments which are much 
lighter than the primary beam are more accurately represented by the approximate solutions 
as seen in figures 5.3 and 5.4. This improvement results from the greater number of terms in 
the summation of equation (3.44). This leads us to believe that the higher order terms in the 
perturbation series can be adequately represented by the approximation in equation (5.8). This 
is especially true because higher order terms in many applications are only small corrections. 

5.2. Approximate Realistic Beams 

Approximate solutions for realistic ion beams may be found by using a super- 
position of the approximate monoenergetic beam solutions. The incident ion beam is taken as 
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where E 0 is the nominal beam energy, and A is related to the half-width at half-maximum. 
The first term is then as before 
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where R M (Em) = Rm(E) + x • 0ne similarly arrives at 
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where 
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and Ej u and E jt are given by equations (5.3) and (5.4). Additional computation yields 
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where Ej u and Eji are given by equations (5.6) and (5.7), and Ei and E u are given in 
equations (4.8) and (4.9). The remaining higher order terms are taken as 


4>f(x,E)= Y, 


a jj n _ j Oj i , M g( j, j„-i , • . . , j \ i M) er f [(E m - 2A] erf [(Ei E u )/\/ 2A 
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where Eji and E ]U are given by equations (5.2) and (5.3) and E u and Ei are given by 
equations (5.13) and (5.14). 

These approximate equations for realistic ion beams are given in figures 5.5 to 5.8 and should 
be compared with the more exact formulae given in section 4. The primary ion beam is taken 
as 20 Ne at 600 MeV /amu with a 0. 2-percent energy spread. The 17 0 flux first collision term is 
shown in figure 5.5 for the two formalisms. The effect of the beam energy spread is seen as a 
rounding of the spectrum at the edges compared with the monoenergetic case in figure 5.1. The 
second collision term is shown in figure 5.6. The approximate second collision term improves for 
the lighter fragments as seen in figures 5.7 and 5.8. Higher order collision terms are expected 
to be more accurate because of the large number of combinations of contributing ion terms. 




Figure 5.1. 17 0 flux spectral term according to first collision term and approximate first collision 
term. 



Figure 5.2. 7 0 flux spectral term according to second collision term and approximate second 
collision term. 
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Figure 5.3. 16 0 flux spectral term according to second collision term and approximate second 
collision term. 



Figure 5.4. 14 N flux spectral term according to second collision term and approximate second 
collision term. 
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III 



E, MeV/amu 


Figure 5.5. flux spectral term for energy spread ^°Ne beam according to first collision term 
and approximate first collision term. 



Figure 5.6. 17 0 flux spectral term for energy spread 20 Ne beam according to second collision 
term and approximate second collision term. 
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Figure 5.7. 16 0 flux spectral term for energy spread 20 Ne beam according to second collision 
term and approximate second collision term. 



Figure 5.8. 14 N flux spectral term for energy spread 20 Ne beam according to second collision 
term and approximate second collision term. 
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6. Recommended Methods 

An energy-independent theory has been used to show that the perturbation expansion up to 
the double collision term is adequate for all fragments whose mass is near that of the projectile. 
This is why the three- term expansion was able to explain the Bragg curve data for 2 ®Ne beams in 
water with reasonable accuracy (ref. 7). As a starting point for the calculation of the transition 
of heavy ion beams in materials, the use of the three-term expansion of sections 3 and 4 can 
be further corrected by use of the approximate higher order terms given in section 5. As an 
example of such a procedure, we give results for 20 Ne beams at 600 MeV/amu in water. The 
results are shown in figures 6.1 to 6.6 as successive partial sums of the perturbation series. The 
solid line is the first collision term. The dashed curve includes the double collision terms. The 
long-dash-short-dash curve includes the triple collision term and can hardly be distinguished 
from the long-dash-double-short-dash curve which includes the quadruple collision terms. The 
results for penetration to 20 cm of water are shown in figures 6.1 to 6.6. The monoenergetic 
beam results for ^O, and * 2 C are given in figures 6.1 to 6.3, respectively. The double 
collision term is seen to be always an important contribution. The triple collision term shows 
some importance for 12 C, while higher order terms are negligible. Similar results are shown in 
figures 6.4 to 6.6 for an energy spread of 0.2 percent. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
May 11, 1989 
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Figure 6.L Sequence of approximations of 17 0 flux spectrum after 20 cm of water for first-order, 
second-order, third-order, and fourth-order theories. 



Figure 6.2. Sequence of approximations of 16 0 flux spectrum after 20 cm of water for first-order 
second-order, third-order, and fourth-order theories. 
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Figure 6.3. Sequence of approximations of 12 C flux spectrum after 20 cm of water for first-order, 
second-order, third-order, and fourth-order theories. 



Figure 6.4. Sequence of approximations for energy spread solution of 17 0 flux spectrum after 
20 cm of water for first-order, second-order, third-order, and fourth-order theories. 
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Figure 6.5. Sequence of approximations for energy spread solution of 16 0 flux spectrum after 
20 cm of water for first-order, second-order, third-order, and fourth-order theories. 



Figure 6.6. Sequence of approximations for energy spread solution of 12 C flux spectrum after 
20 cm of water for first-order, second-order, third-order, and fourth-order theories. 
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